The pseudogap Fermi-Bose Kondo model 
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We consider a magnetic impurity coupled to both fermionic quasiparticles with a pseudogap 
density of states and bosonic spin fluctuations. Using renormalization group and large-7V calculations 
we investigate the phase diagram of the resulting Fermi-Bose Kondo model. We show that the 
Kondo temperature is strongly reduced by low-energy spin fluctuations, and make connections to 
experiments in cuprate superconductors. Furthermore we derive an exact exponent for the critical 
behavior of the conduction electron T matrix, and propose our findings to be relevani for certain 
scenarios of local quantum criticality in heavy-fermion metals. 



The interplay of quasiparticles and coUective low- 
energy excitations is a central theme in the physics of 
strongly correlated materials. Impurities have proven to 
be a powerful probe for investigating the bulk behav- 
ior of such Systems. A paradigmatic model describing 
the interaction of impurity degrees of freedom with both 
quasiparticles and coUective spin fluctuations is the so- 
called Fermi-Bose Kondo model P, 0, consisting of a 
local impurity spin, S, coupled to spin- i fermions and 
spin-1 vector bosons. 

The Fermi-Bose Kondo model has recently been ana- 
lyzed for the case of a metallic fermion density of states 
(DOS), and a gapless bosonic spectrum representing 
magnetic order parameter fluctuations at a bulk quantum 
critical point in d dimensions P, 1^ ■ For d < 3 the model 
shows a boundary quantum phase transition between a 
Kondo-screened phase and a bosonic fluctuating phase 
with universal local spin correlations. Both the critical 
and the bosonic fluctuating fixed points are characterized 
by singular thermodynamic properties, connected with 
local non-Fermi liquid behavior. 

The purpose of this Letter is to generalize this analysis 
to a fermionic bath with a pseudogap DOS, Pc(^) c>c 
and to discuss Kondo screening for the case of a flnite 
spin gap. As, in the bosonic bath. Using renormaliza- 
tion group (RG) methods, we shall establish the phase 
diagram of the pseudogap Fermi-Bose Kondo model, and 
furthermore derive an exact relation between the anoma- 
lous exponent of the fermion T matrix and the expo- 
nent of fermionic bath DOS, r, which holds at any fixed 
point with a flnite coupling between fermions and impu- 
rity spin. We shall apply our model to impurity momcnts 
in cuprate superconductors, where the fermionic DOS 
obeys r = 1, and show that Kondo screening is strongly 
suppressed by coUective spin fluctuations, in agrccmcnt 
with NMR experiments. 

Recently, the Fermi-Bose Kondo model has been pro- 
posed m to describe a "local quantum phase transition" 
in alloys Hke CeCu6-xAux This modeUing is based 
on an extended dynamical mean-field theory where 
a lattice model is mapped onto a self-consistent single- 
impurity model. Based on our T matrix analysis, we 
will argue that no anomalous power law can occur in the 



fermion spectrum within such a scenario of local critical- 
ity. 

Model. The Fermi-Bose Kondo Hamiltonian for a spin- 
i impurity can be written a,s H — Hi + H\y + i?imp, with 



-ffimp = JkS ■ Sb + 70^ • 00 • 



(1) 



Here, S denotes the impurity spin operator, and Hi de- 
scribes a system of free fermions, Hf = ^Z^ko- ^kcj^-^Cko- 
in Standard notation, so = X] 



kk' 



„t 



f?crcr'Ck'cr' IS thc COU- 



duction band spin operator at the impurity site ro = 0. 
The fermionic density of states foUows a power law at low 
energies, Pc(^^) = X]k '^('^ ~ ^k) oc iVo|u;|'', which includes 
the cases of a metal (r = 0) and a d-wave superconductor 
(r = 1). In thc simplest case, the bosonic bath consists 
of free vector bosons, H\, — X]qQ '^q^qQ^qai describing 
coUective spin fluctuations in the host material, with a 
dispersion cj^ = m? + c^q^. At zero temperature, the 
mass m coincides with the bulk spin gap As. A gap 
As > O describes a quantum paramagnet. As = O cor- 
responds to a bulk quantum critical point between the 
paramagnet and an antiferromagnetically ordered phase; 
the momentum q is measured relative to the ordering 
wavevector Q. The field 0o in Eq- lQl represents the lo- 
cal orientation of the antiferromagnctic order parameter, 
it is given by 0oa = Z]q(^qa + ^-qa)/ V^q/>^i whcrc J is 
the bulk exchange constant. It is useful to deflne a spec- 
tral density, P(j,{uj), of the Bose field (p. At T = O and for 
As = O, it foUows p^{uj) cx K§ sgn(w)|cj|i-% = J/d^, 
in (i = 3— e dimensions. We note that in dimensions d < 3, 
boson self-interactions are relevani in the RG sense (in 
the absence of Landau damping), and we will comment 
on their effect below. Our analysis will be restricted to 
thc SU(2)-invariant model; the effect of spin anisotropy 
for r = has been discussed in Ref. 2. 

In the absence of the bosonic bath, 7 = 0, and for r = 
the model (Q reduces to the well-known Kondo model |^ , 
which cxhibits screening of thc magnetic momcnt below a 
characteristic (Kondo) temperature, Tr-. The pseudogap 
generalization, r > O, shows a zero-temperature phase 
transition between a free local moment and a screened 
moment as the Kondo coupling is varied 0, S 0] ■ On 
the othcr hand, in the absence of the fermionic bath. 





FIG. 1: Schematic RG flow diagram for the pseudogap Fermi- 
Bose Kondo model for small e and r, and As = 0. The axes 
denotes the couphngs to the bosonic (7) and fermionic (j) 
baths. The model has three (meta)stable fixed points (solid 
dots): the local moment (LM), bosonic fluctuating (B-FL), 
and fermionic strong-coupling (F-SC) fixed points; the open 
dots are the critical fixed points F-CR and FB-CR; for details 
see text. The thick line corresponds to a line of continuous 
boundary phase transitions. In the metallic case, r — O, the 
LM and F-CR fixed points merge |1I3. 

J K — O, the Bose-Kondo model deseribes an impurity 
moment in an insulating quantum antiferromagnet [loj . 

RG analysis. At the bulk quantum critical point we 
have Ag = O, and both fermionic and bosonic baths obey 
a power law DOS. The impurity behavior can be ana- 
lyzed using field-theoretic RG methods, together with an 
expansion in e and r, similar to Refs.IlLIMIlCll 

We employ a pseudo-fermion representation of the im- 
purity spin i, S = flSaa' fa- FoUowing the Standard 
procedure we introduce renormalized fields = 

■y/ Zffna-, (j) a = ^/Z4)Rcn and dimensionless couplings 

i^o7o = /i'/'[^7/(^/V^)]7, NqJk = ti'^{Zj/Zf)j, 
where /i is a renormalization energy scale. The deriva- 
tion of the RG equations has been presented in Ref. for 
non-intcracting (j) bosons (Z = 1 in this case). A similar 
procedure for the present case with a fermionic pseudo- 
gap DOS, using the minimal subtraction scheme, gives 
the foUowing result for the RG beta functions: 

/3(7) = -y+7 -7 + — , 

m = n-/ + y +J7'-J7'- (2) 

The fuU two-loop result for the case of interacting bosons 
will be given elsewhere. The equations (01 reduce to the 
cases of the pseudogap Kondo model for 7 = O 6] and the 
metallic Bose-Fermi Kondo model for r = O j^. Besides 
the decoupled (local moment - LM) fixed point, 7* = j* = 
O, the above flow equations have three non-trivial fixed 
points which are perturbatively accessible: 
B-FL: 7* = I -I- j* = O corresponds to a stable 
fixed point where the fermionic Kondo coupling is irrel- 
evant, and the coupling between impurity and bosonic 
bath leads to universal local- moment fluctuations (lO] : 



FIG. 2: Feynman diagrams for the Green's function Gt en- 
tering the conduction electron T matrix. FuU/dashed/wiggly 
lines denote j jcjc^ propagators. The double line is the full 
/ propagator, the shaded area the full vertex, and the open 
dots are the external sources. 
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F-CR: j* = r -f 7* = O is an unstable fixed point 
controUing the transition between the decoupled (free) 
moment and the Kondo-screened moment phases; and 
FB-CR: 7*2 = I + ^ _ ^ _ j* = r + I - this is an 
unstable fixed point between the bosonic fluctuating and 
the Kondo-screened phases. 

For large enough j, 7 is irrelevant, and the model re- 
duces to the fermionic pseudogap Kondo model which 
shows strong-coupling phases with screening. Numerical 
work has shown that, depending on the value of the 
exponent r and the presence or absence of particle-hole 
symmetry, two different strong-coupling fixcd points ex- 
ist. As these are not accessible within the present pertur- 
bation expansion, we will not distinguish between them, 
and simply denote the fermionic strong-coupling fixed 
point, reached for large j and small 7, with F-SC. The 
flow diagram resulting from this discussion is shown in 
Fig. ni In the presence of bosonic interactions the RG 
flow is similar, and the location of the fixed points is 
modified only at order e^. We expect the general struc- 
ture of the phase diagram, Fig. ^ to be valid for all e 
and r values between O and 1 - this is e. g. the result of 
large- iV calculations (see below); for r > i particle-hole 
asymmetry is needed to reach the F-SC fixed point. Fi- 
nally, a finite host spin gap, > O, will cut-off the RG 
flow of 7 at the scale As, but there is still a transition as 
function of Jk controlled by the F-CR fixed point. 

T matriz. An important quantity is the conduction 
electron T matrix, describing the scattering of the c elec- 
trons off the impurity. It is useful to define a propagator, 
Gt, of the composite operator — /J/o-'Ccr', such that 
the T matrix is given by T(lo) — J]^Gt{(^)- The dia- 
grams for Gt in perturbation theory are shown in Fig. |2] 
At the RG fixed points, we expect a power law behavior of 
the T matrix spectral density, with an anomalous expo- 
nent rjT- (Note that at tree level Gt oc w'" in the present 
pseudogap problem.) Within the Standard scheme 
we can obtain the Gt renormalization factor Zt in an 
expansion in 7 and j by evaluating the lowest-order dia- 
grams and demanding the resulting expression to be free 
of poles at energy /j,. From Zt we obtain a perturbative 
expression for rjT according to rjT = d{\ii Z t)/ d (In 11), 
evaluated at j = j* and 7 = 7*. 
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However, it is possible to obtain an ea;ac<result for rjT, 
valid to all orders in perturbation theory. The argument 
parallels the one for anomafous exponent of the focal sus- 
ceptibility in Refs. '5,'lo'|: Inspecting the diagrams in the 
upper panel of Fig.|21shows that all are obtained from the 
first diagram with full j vertices and full f electron prop- 
agators, i. e., there are no additional processes involving, 
e. g., the conduction electron line. This gives the relation 
Zy"^ — Zj/Zj Using the definition of Zj and taking 
the derivative w. r. t. (In/i) at fixed bare Jk, we obtain 
O = —r + (3{j)/j — Vt/'^- Thus, at any RG fixed point 
with finite, non-zero j*, we find the exact result 

This result applies to both critical fixed points (F-CR and 
FB-CR), and relies only on the fact that the c electrons 
do not have a self-interaction being relevant in the RG 
sense, i. e., it holds for interacting ^ bosons as well. 

Remarkably, the critical behavior Tiyi) cx w"'' has 
been numerically found in studics of the pseudogap 
Kondo model, i. e., at the F-CR fixed point 0, 0|- It 
is also obtained within a dynamic large-A'^ analysis of the 
same problem 

Application to cuprate superconductors. In supercon- 
ducting cuprates it is reasonable to separate the low- 
energy degrees of freedom into Bogoliubov quasiparti- 
cles around the nodes of the d-wave superconducting or- 
der parameter, and antiferromagnetic spin fluctuations 
around wavevector Q — (7r,7r), as manifest in the reso- 
nance peak in neutron scattering ; momentum conser- 
vation prohibits a low-energy coupling between fermionic 
and bosonic degrees of freedom. A phenomenological de- 
scription of cuprates in terms of a BCS superconduc- 
tor plus low-energy coUective spin excitations [T^ has 
been successful in describing a number of recent neutron 
scattering and STM experiments. Therefore, impurity 
moments in cuprate d-wave superconductors can be ex- 
pected to be well described by the pseudogap Fermi-Bose 
Kondo model, with fermionic DOS exponent r — l. 

Experimentally, signatures of Kondo physics have been 
observed in cuprates doped e. g. with Zn replacing Cu 
[ifil Ir^ - here the non-magnetic Zn induces a quasi-free 
magnetic moment in its vicinity, which can be related to 
the properties of the parent Mott insulator . In par- 
ticular, NMR experiments indicate Kondo screening of 
Zn-induced moments in optimally doped YBa2Cu307_5 
below the superconducting Tc, but show that screening 
is suppressed in underdoped samples far, theo- 

retical analysis has been restricted to either a moment 
interacting with bulk spin fluctuations |lO| or a moment 
interacting with Bogoliubov quasiparticles 0, lisl ligj . 

Here we propose that the increasing strength of anti- 
ferromagnetic fluctuations is the most important contri- 
bution to the Tk suppression in underdoped cuprates. 
To this end, we calculate the Kondo temperature in the 




Spin gap A, [meV] 

FIG. 3: Characteristic energy scale T* as function of the 
bosonic spin gap As, calculated for a fermionic host cor- 
responding to a cuprate superconductor p^ . using large-A'' 
techniques. Solid: Jk =0.27 eV, U = -0.15 eV, 70 = 0.15 
eV; Dash-dot: Jk =0.38 eV, (7 = 0, where U is the additional 
impurity scattering potential. In comparison with NMR ex- 
periments. T* is equivalent to the Kondo temperature below 
which the impurity moment is screened; the horizontal axis 
can roughly be identified with doping (up to As ~ 50 meV) 
fl^ . Note that the change in the fermionic band upon doping 
has little influence on T* ; an increase in the superconducting 
gap, Ao, with underdoping will further suppress T*. 

pseudogap Fermi-Bose Kondo model for a realistic band 
structure in Hf and a spin fluctuation spectrum in 
with varying spin gap A^. In modelling the experiment, 
the important ingredient is that Ag decreases with un- 
derdoping, as observed in neutron scattering [l^ - from 
the RG analysis it is clear that decreasing Ag can drive 
the impurity from the F-SC to the LM phase. For Hf we 
use a Standard d-w&ve BCS model with a gap function 
Ak = (Ao/2)(cos fca; — cosA:^). 

To obtain quantitative results, we employ a large-A^ 
approach were the spin symmetry is generalized from 
SU(2) to S UfiV)- T he methodology is similar to earlier 
work [5I Hol ll^ Hoj. details will be published elsewhere. 
While this method has artifacts near the quantum- 
critical points discussed above, and is not quantitatively 
accurate, we expect it to capture the qualitative low-T 
physics. The T = O results allow to extract a low-energy 
scale, namely the position of the peak in the T matrix, 
ksT*, which can be identified with the Kondo temper- 
ature up to a numerical factor of order unity jl3l Il9j . 
Fig. 13 shows sample results for this quantity as function 
of the bulk spin gap, Ag. 

These results indicate that antiferromagnetic fluctua- 
tions with a small spin gap A^ are effective in suppress- 
ing Kondo screening. We note that an increasing su- 
perconducting gap, Ag, upon underdoping can have a 
similar effect in the pseudogap Kondo model How- 
ever, this mechanism is less likely to be universal for all 
cuprates, as a strong doping dependence of Ao has been 
observed in Bi2Sr2CaCu208+5 but not in YBa2Cu307_5. 
In contrast, spin fluctuations grow upon underdoping in 
all cuprates. Therefore, we propose the suppression of 
Kondo screening as observed in NMR to be well ex- 
plained within the pseudogap Fermi-Bose Kondo model 
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with a doping-dependent host spin gap. 

Application to heavy-fermion metals. Motivated by 
neutron scattering experiments on the heavy-fermion 
compound CeCug_xAux 0, which indicate momentum- 
independent critical dynamics at an antiferromagnetic 
ordering transition, Si et al. proposed a self-consistent 
version of the Fermi-Bose Kondo model as a candidate 
to deseribe such critical behavior within an extended dy- 
namical mean-field approach (EDMFT) 0]. In this sce- 
nario, the critical point of the lattice model is mapped 
onto the critical point (BF-CR) of the impurity model. 

The self-consistency within EDMFT has to be achieved 
for both fermionic and bosonic baths. So far, most dis- 
cussions assumed r = O and focussed primarily on the 
bosonic part 0, 0] ■ Here we want to draw attention to 
the fermionic part as well. Our analysis above shows 
that a fermionic bath with \uj\'' spectrum leads to a im- 
purity T matrix with power-law divergence \uj ® at 
the critical point. Now the question is whether DMFT 
self-consistency in the fermionic bath of a Kondo or An- 
derson lattice model can be achieved with such a power 
law. Consider the DMFT self-consistency relation for an 
Anderson lattice model, where the T matrix corresponds 
to F2Gi,„p ^5]: 



G'°^z) = J de— 



Po(e) 



e/ -S/ (z) 
-A(z)-I]/(z)] 



(4) 



^imp (^) 



in Standard notation, with /9Q(e) the DOS of the bare 
conduction band of the lattice model, S-^ the intcracting 
impurity self-energy, and A (z) the effective hybridiza- 
tion function for the impurity model. Now, Gimp oc uj^^ 
implies Im[A(w) -|- iuj)] oc — As Yi^ is identi- 
fied with the lattice self-energy within DMFT, we require 
ImE-^ < O, and both ImA(a;) and ItuTj^ {uj) have to van- 
ish as — Iw^- Using the ansatz Y^ {z) — —Sf + Az^ and 
expanding /o3(e) around e = O, the integral in the first 
line of ^ can be performed in the limit of small z, and, 
surprisingly, the coefficient of the leading z^^ term in 
G^°'^{z) vanishes. 

Therefore, critical-point self-consistency in the 
fermionic sector of EDMFT for an Kondo or Anderson 
model cannot be achieved with an exponent r > 0. If 
a local critical point within EDMFT exists, it is thus 
characterized by r = O, and in this case logarithmic 
corrections in the T matrix are expected, leading to a 
logarithmic temperature dependence of the resistivity. 
In addition, we note that any such critical point will 
be characterized by an extensive ground state entropy 



(this holds for arbitrary r), which should have strong 
signatures in specific heat experiments. 

Conclusions. We have analyzed a Fermi-Bose Kondo 
model, describing an impurity spin coupled to both 
fermionic quasiparticles with a pseudogap DOS and col- 
lective mode bosons, using renormalization group and 
large-A^ techniques. We have argued our findings to be 
relevant for cuprate superconductors as well as for sce- 
narios of local criticality in heavy-fermion metals. 
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